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Abstract 

We perform a gfobal analysis of the latest solar neutrino data including the Solar 
Neutrino Observatory (SNO) result on the CC— event rate. This result further favors 
the LMA MSW solution of the solar neutrino problem. The best fit values of parame- 
ters we find are: Am^ = (4.8 - 5.0) • IQ-^ eV^, tan^ d = 0.35 - 0.38, fs = 1.08 - 1.12, 
and fhep = 1 — 4, where /b and fhep are the boron and hep— neutrino fluxes in units 
of the corresponding fluxes in the Standard Solar Model (SSM). With respect to this 
best fit point the LOW MSW solution is accepted at 90 % C.L.. The Vacuum oscilla- 
tion solution (VAC) with Am^ = 1.4 • 10"^*^ eV^, gives good fit of the data provided 
that the boron neutrino flux is substantially smaller than the SSM flux (fs ~ 0.5). 
The SMA solution is accepted at about 3a" level. We find that vacuum oscillations 
to sterile neutrino, VAC(sterile), with fs ~ 0.5 also give rather good global fit of the 
data. All other sterile neutrino solutions are strongly disfavored. We check the quality 
of the fit by constructing the pu/Z-Oj^ diagrams of observables for the global solutions. 
Maximal mixing is allowed at 3a level in the LMA region and at 95 % C.L. in the 
LOW region. Predictions for the day-night asymmetry, spectrum distortion and ratio 
of the neutral to charged curent event rates, [NC]/[CC], at SNO are calculated. In 
the best fit points of the global solutions we find: « (7 - 8)% for LMA, ~ 3% 

for LOW, and (2 — 3)% for SMA. In the LMA region the asymmetry can be as large 
as 15 - 20 %. Observation of ^dg^ > 5% wih further favor the LMA solution. It wih 
be difficult to see the distortion of the spectrum expected for LMA as well as LOW 
solutions. However, future SNO spectral data can significantly affect the VAC and 
SMA solutions. We present expected values of the BOREXINO event rate for global 
solutions. 



Pacs numbers: 14.60.Lm 14.60.Pq 95.85.Ry 26.65.+t 



1 Introduction 



With the SNO result [T|, H on the CC— rate, we now, for the first time ever, have more than 
3(7 evidence of flavor conversion of solar neutrinos. "Smoking guns" have indeed started 
to smoke. The statement, made first in the original publication by the SNO collaboration, 
is based on difference of the boron neutrino fluxes determined from the charged current 
{CC—) event rate in the SNO detector, and the z/e— scattering event rate obtained by 
the SuperKamiokande (SK) collaboration (and confirmed, albeit with smaller statistical 
significance, by SNO). Somewhat paradoxically, the absence of significant distortion of the 
boron neutrino spectrum at SuperKamiokande adds to the strength of this conclusion. 

In general, the difference between the signals in the SNO and SK detectors can be due 

to: 

1) appearance of the i/^, Vt-— flux, or/and 

2) distortion of the neutrino energy spectrum. If the suppression of the boron neutrino 
flux (due to some neutrino transformations) increases with neutrino energy, then the higher 
sensitivity to higher energies of the CC— reaction in SNO as compared with v — e scattering 
in SK would explain the lower rate in SNO. 

In fact, both reasons imply neutrino conversion. However, an absence of strong distortion 
of the boron neutrino spectrum, as found by SuperKamiokande and independently confirmed 
by SNO, leads to the conclusion that the main reason for the difference is the appearance 
of the i/^, z/,-— flux. 

From the SNO result and its comparison with the SK data one can immediately draw 
several conclusions: 

• Strong deficit of the z/g— flux with respect to the Standard Solar Model (SSM) pre- 
dictions is found. 

• There is a strong evidence that Ve from the sun are converted to v^iVt- 

• There is no astrophysical solution of the solar neutrino problem. (For recent quanti- 
tative analysis see @]). 

• Solutions of the solar neutrino problem based on pure active - sterile conversion, 

i-'si are strongly disfavored. 

• More than half of the original z/g— flux is transformed to neutrinos of different type z/^, 
Vri and possibly Vg-, that is, the z/g— survival probability 

P < 1/2. (1) 

In fact, in assumption of pure active transition one gets P = 0.334 ± 0.22 which is 
just ~ Icr below 1/2. Clearly, if the sterile neutrino component is present in the flux. 
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the survival probability is even smaller. The inequality ([l|), if confirmed, will have 
crucial implications for further experimental developments in the field, as well as for 
fundamental theory. 

We know now with a high confidence that electron neutrinos produced in the center 
of the Sun undergo flavor conversion. However the specific mechanism of the conversion 
has not yet been identified. As we will see, only some extreme possibilities are excluded 
by adding the SNO result to the previously available solar neutrino data. A number of 
solutions still exist. 

The SNO result changes the status of specific solutions to the solar neutrino problem. 
The changes can be immediately seen by comparing predictions for the CC— event rate 
from global solutions found in the pre-SNO analysis 0, 0| with the SNO result. There are 
four solutions, three active and one sterile, which give predictions close to the SNO result: 

^Sc''^ IIS = 0.347 ±0.029. (2) 

SSM ^ ' 

1) . The Large Mixing Angle MSW solution (LMA): the 3ct predictions interval, Rcc = 
0.20 - 0.41, covers the SNO result {Rcc = [CC] from |]). The best fit point from the 
pre-SNO analysis, Rcc = 0.31, is slightly (~ Icr) below the central value given by SNO. 
Therefore, SNO shifts the region of the LMA solution and the best fit point to larger values 
of mixing angles which correspond to larger survival probability. 

2) . The low Am^ (LOW) solution: the interval of predictions, Rcc = 0.36 — 0.42, is 
above the SNO result. In the best fit point Rcc is 2(j (experimental) higher the central 
SNO value. Therefore this solution is somewhat less favored, and SNO tends to shift the 
allowed region to smaller values of 6 which correspond to smaller survival probability. 

3) . The vacuum oscillation solutions (VAC): the expected interval, Rcc = 0.33 — 0.42, 
is also covering the SNO range, and the best fit value, Rcc = 0.38, is just la above the 
central SNO result. Therefore, SNO improves the status of this solution. 

4) . Vacuum oscillations to sterile neutrinos, VAC(sterile): the predicted interval, Rcc = 
0.36 — 0.41, with the best fit point 0.39, is only slightly higher than for the VAC(active) 
case. Here low rate predicted for SNO is due to difference of the thresholds in SNO and SK 
and to the steep increase of the suppression with neutrino energy. This solution survives 
the first SNO result. 

Other solutions look less favorable in the light of the SNO result. In particular, the 
small mixing angle MSW (SMA) solution interval, Rcc = 0.37 — 0.50, (especially its best fit 
point, 0.46, from the pre-SNO analysis) is substantially above the SNO rate. This solution 
is further disfavored by SNO. In fact, SMA predicts opposite energy dependence of the 
suppression on the neutrino energy to the one SNO prefers. 

The Just-so^ (active), SMA (sterile) and Just-so^ (sterile) solutions predict practically the 
same rate for the SK and SNO and therefore are strongly disfavored. 
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The statements above are supported by detailed quantitative analysis of the solar neu- 
trino data we have performed in this paper. The results of this analysis are described in 
the rest of this paper. Previously some of the same conclusions have been reported in 
0, |, |, |T3; see also related studies [[11], |12|, [13], 0. 

The paper is organized as follows: In Section 2 we present the results of the global anal- 
ysis of the solar neutrino data. We construct the pull-off diagrams of available observables 
for the found global solutions. In Section 3 we consider specifics of global solutions and 
their implications. We evaluate the quality of the data fit in each of the global solutions. 
We calculate predictions for the day-night asymmetry and spectrum distortion in Section 4. 
Finally, prospects for identifying the solution to the solar neutrino problem are discussed in 
Section 5. 



2 Global analysis and pull-ofF diagrams 

We describe here the results of the global analysis of the solar neutrino data. We find global 
solutions (sets of oscillation parameters and solar neutrino fluxes which explain the solar 
neutrino data) and determine the goodness of fit (g.o.f.) in the best fit points. We perform 
diagnostics of the global solutions by checking their stability with respect to variations in 
the analysis and to uncertainties in the original solar neutrino fluxes. To check the quality 
of the fit we construct the pull-off diagrams for the observables. 



2.1 Features of the analysis 



We follow the procedure of the analysis developed in our previous publications |T^, [T^, 0, Q 
in collaboration with John Bahcall. We also describe some additional features, which should 
be taken into account when comparing our results with those obtained by other groups. 

In our global analysis we use: (i) the Ar— production rate from the Homestake experi- 
ment (ii) the Ge— production rate from SAGE |jl9|, (iii) the combined Ge— production 



rate from GALLEX and GNO 0, 0, (iv) the CC-event rate measured by SNO 0, (v) 
the Day and Night energy spectra measured by SuperKamiokande 0. 

Following the procedure outlined in p[ we do not include the total rate of events in 
the SK detector, which is not independent from the spectral data. In our standard global 
analysis, we use 4 rates and 38 spectral data, a total of 42 data points. The number of free 
parameters and the number of d.o.f. are different in different analyses and we will specify 
them later. 



The solar neutrino fluxes are taken according to SSM BP2000 |22| with the corrected 



(due to improved measurement of the solar luminosity) boron neutrino flux Ff'^^^ = 5.05-10^ 
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cm ^ c ^. We denote by Jb and fhep the fluxes of the boron and /lep— neutrinos measured 
in units of the BP2000 fluxes. 

The analysis of the data is performed in terms of two neutrino mixing characterized by 
the mass squared difference Am^ and the mixing parameter tan^ 6. We consider conversion 
into pure active or pure sterile neutrinos. 

We perform the test of various oscillation solutions by calculating 



^global ^rate ' "X. spectrumi 

where Xrate 'ydspectrum ^'^^ coutributious from the total rates and from the Su- 
per Kamiokande day and night spectra correspondingly. Each of the entries in Eq.(^ is 
a function of the four parameters (Am^, tan^ 6, fs and fhep)- It is an important feature of 
our approach that at each step of the minimization process these parameters are kept 
the same in each of the two entries on the right hand side of Eq. (j^) . 



2.2 Free flux fit 



We perform the flt to the experimental data treating the boron neutrino flux, fs, and the 
/lep— neutrino flux, f^ep, as free parameters. There are several reasons to consider Jb and 
fhep as free parameters (for earlier work in this context see [^ ): 

1. These fluxes have the largest uncertainties in the SSM (see however, ||2^ ). 

2. The goal of the solar neutrino studies is to flnd directly from the solar neutrino data 
both the oscillation parameters and the neutrino fluxes. The free flux analysis is the way 
to achieve this goal. 

3. Comparing the fluxes fhep and found from the free flux analysis with the SSM 
fluxes we can estimate the plausibility of the fit and the reliability of the solutions. Clearly, 
strong deviation of the fluxes for a given solution from the SSM values will indicate certain 
problems with either the solution or the SSM predictions. 

4. Last, but not least, this flt is the most conservative one regarding the exclusion of 
certain scenarios. 

Thus, together with Am^ and tan^ 6, there are 4 free parameters and therefore 42 - 4 = 
38 d.o.f. in the x^'ht. In Table |I] we show the best flt values of the parameters Am^, tan^ 6, 
fs, fhep for different solutions of the solar neutrino problem. We also give the corresponding 
values of xLin goodness of the flt. Several remarks are in order. 

The absolute minimum: = 29.2 is in the LMA region. Such a low for 38 d.o.f. 
is mainly due to the small spread of the experimental points in both the day and the night 
SK spectra. Similar situation is realized for the LOW solution. The SuperKamiokande day 
and night spectra can be especially well described by solutions which predict a bump in the 
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Table 1: Best-fit values of the parameters Am^, tan^ ^, fs and fhep from the free fiux 
analysis. The minimum cind the corresponding g.o.f. are given in the last two columns. 
The number of degrees of freedom is 38: 4 rates (Homestake, SAGE, Gallex/GNO, SNO) 
+ 38 SK spectra points - 4 parameters. 



Solution 




/eV 


tan2 6' 


fB 


fhep 


Anwn 


g.o.f. 


LMA 


4.8 


X 


10- 


-5 


0.35 


1.12 


4 


29.2 


0.85 


VAC 


1.4 


X 


10" 


-10 


0.40 (2.5) 


0.53 


6 


32 


0.74 


LOW 


1.1 


X 


10" 


-7 


0.66 


0.88 


2 


34.3 


0.64 


SMA 


6.0 


X 


10 


-6 


1.9 X 10-3 


1.12 


4 


40.9 


0.34 


Just So^ 


5.5 


X 


10- 


12 


1.0 


0.44 





45.8 


0.18 


VAC (sterile) 
Just So^ (sterile) 


1.4 


X 


10- 


-10 


0.38 (2.6) 


0.54 


9 


35.1 


0.60 


5.5 


X 


10" 


-12 


1.0 


0.44 





46.2 


0.17 


SMA (sterile) 


3.8 


X 


10" 


-6 


4.2 X 10-^ 


0.52 


0.2 


48.2 


0.12 


LMA(sterile) 


1.0 


X 


10^ 


-4 


0.33 


1.14 





49.0 


0.11 


LOW(sterile) 


2.0 


X 


10- 


-8 


1.05 


0.83 





49.2 


0.11 



survival probability aX E = Q — 8 MeV and a dip ati?=10 — 11 MeV. This is the case 
of VAC solutions (both active and sterile) with large /lep— fiux. It is for this reason VAC 
solutions have high goodness of the global fit. According to the Table || the VAC (sterile) 
solution is in the third position after LMA and VAC (active) and its fit is even better than 
the one of LOW solution. The VAC(sterile) solution has, however, a number of problems 
which we will discuss in the next section. 

For LMA, LOW and SMA solutions values of J'b, agree with the SSM predictions within 
la theoretical uncertainty (~ 18 %). All VAC solutions and SMA(sterile), as well as Just-so^ 
solutions, appear with a boron fiux which is 3a below the SSM boron neutrino fiux. For the 
/lep— neutrino fiux the VAC and SMA solutions imply significant (factors 4-6) deviation 
from the central SSM value. The VAC(sterile) requires as large as 9 SSM /lep— neutrino 
fiux. 

For VAC solutions the matter effect is negligible and in the Table 1, two "symmetric" 
values of tan^ 6: tan^ di = 1/ tan^ 62 correspond to mixing in normal and dark (in brackets) 
sides of the parameter space. 

Even within this conservative analysis the LMA(sterile) and LOW(sterile) solutions give 
a very bad fit and in what follows we will not discuss them. 

Just-so^ (active) and Just-so^ (sterile) give very similar description of the data. So, we 
will present results for one of them. 

In Fig. H we show contours of constant (90, 95, 99, 99.73 %) confidence level with respect 
of the absolute minimum in the LMA region. Following the same procedure as in the 
contours have been constructed in the following way. For each point in the Am^, tan^ 6 
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Table 2: Best-fit values of the parameters Am^, tan^ 9 and fs from the global analysis with 
fhep = 1- The minimum and the corresponding g.o.f. are given in the last two columns. 
The number of degrees of freedom is 39: 4 rates (Homestake, SAGE, Gallex/GNO, SNO) 
+ 38 SK spectral points - 3 parameters. 



Solution 






tan2 6' 


fB 


A mm 


g.o.f. 


LMA 


5.0 


X 




0.36 


1.1 


30.1 


0.85 


VAC 


1.4 


X 


10-10 


0.363 (2.7) 


0.54 


33.4 


0.72 


LOW 


1.1 


X 


10-7 


0.69 


0.86 


34.5 


0.67 


SMA 


5.5 


X 


10-^ 


1.9 X 10-3 


1.04 


42.2 


0.33 


Just So^ 


5.5 


X 


10-12 


1.0 


0.44 


46.4 


0.19 


VAC (sterile) 


1.4 


X 


10-10 


0.35 (2.9) 


0.55 


36.9 


0.57 


SMA (sterile) 


3.8 


X 


10-*^ 


4.2 X 10-^ 


0.52 


48.6 


0.14 



plane we find minimal value Xmm(^"^^) ^^^^ ^) varying and fhep- We define the contours 
of constant confidence level by the condition 

xLni^m', tan^ 6) = xL^LMA) + Ax^ , (4) 

where x^^^^{LM A) = 29.3 is the absolute minimum in the LMA region and Ax^ is taken 
for two degrees of freedom. 

To clarify the role of the /lep— neutrino flux we present in Table ^ the result of the 
analysis when fhep = 1 and treat only fs as a free parameter. As follows from the Table ^ 
fixing the flux fhep leads to rather small change of the oscillation parameters. At the same 
time, the constraint fhep = 1 lowers a goodness of the fit of solutions which imply large 
value of fhep in the free flux analysis. We get Ax^ = i.g for VAC(sterile), Ax^ = 1.4 for 
SMA and Ax^ = 1.4 for VAC(active). Notice that now VAC(sterile) is shifted to the fourth 
position. 

According to Ref. ?? the calculated hep— neutrino flux has about 20% uncertainty. 
Therefore solutions which require large fhep (2 - 9) are disfavored and the results of the fit 
with fhep = 1 look more relevant. 



2.3 SSM restricted global fit 

In order to check the significance of the SSM restriction on the boron neutrino flux we have 
performed the fit to the data adding to the x^ sum in Eq. (^) the term 
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Table 3: Best-fit values of the parameters Am^, tan^ 6, and fhep from the SSM restricted 
global analysis, fhep is considered as a free parameter. The minimum and the corre- 
sponding g.o.f. are given in the last two columns. The number of degrees of freedom is 
40: 4 rates (Homestake, SAGE, Gallex/GNO, SNO) + 38 SK spectral points + 1 (/s) - 3 
parameters. 



Solution 




/eV^ 


tan2 


Ib 


fhep 


• 

Arrtin 


g.o.f. 


LMA 


5.0 


X 


10" 


-5 


0.36 


1.10 


4.0 


29.5 


0.84 


SMA 


5.5 


X 


10- 


-6 


1.9 X 10-3 


1.10 


4.0 


41.2 


0.33 


LOW 


1.1 


X 


10- 


-7 


0.66 


0.88 


2.0 


34.7 


0.62 


VAC 


4.8 


X 


10- 


10 


1.9 


0.72 


0.0 


35.4 


0.59 


Just So^ 


5.5 


X 


10- 


12 


1.4 


0.44 


1.0 


55.4 


0.03 


VAC (sterile) 


1.4 


X 


10- 


10 


2.63 


0.55 


8.0 


41.3 


0.33 


SMA (sterile) 


4.0 


X 


10- 


-6 


4.8 X 10-^ 


0.54 


2.0 


54.5 


0.04 



where as = 0.18 is the average (of the upper and lower) la theoretical error of the fiux in 
BP2000 model p2|. With the term (|^) included in our global we, in a way, treat the SSM 



prediction for the boron neutrino fiux as independent "measurement" and consider it as an 
additional degree of freedom. This procedure makes sense because a significant contribution 
to the error in SSM determination of the boron neutrino fiux comes from the measurements 
of the p — Be cross-section. 

One "technical" remark is in order. Our approach differs from analyses where the boron 
neutrino flux is taken as a theoretical prediction from the SSM ||^, |^. In the latter case the 
term (^) is absent, the central SSM value is used in the predictions of observables and the 
theoretical errors on the boron neutrino flux are added to the experimental errors. 

In Table ^ and Fig. ^ we present the results of the flt with SSM constrained boron 
neutrino flux. As expected, the most signiflcant changes in comparison with the free flux 
analysis (Table |T]) appear in those solutions and regions of the oscillation parameters which 
imply strong deviation of fs from 1. Our results show that mostly the best flt points, as 
well as the Xmm ^^^1 goodness of the VAC (active). Just-so^, VAC (sterile) and SMA (sterile) 
solutions are affected. Thus, for the VAC(active) solution Ax^j„ = 3.4, for VAC(sterile): 
^xlnin = 6.2. Moreover, the best flt VAC solution shifts to another point of oscillation 
parameters in agreement with results of other groups. There is no signiflcant change of the 
LMA, LOW and SMA best flt points and goodness of the flt. 

In Fig. 1^ we show the contours of constant (90, 95, 99, 99.73 %) confldence level for two 
degrees of freedom with respect to the absolute minimum in the LMA region. 
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Table 4: Best-fit values of the parameters Am^ and tan^ 6 from the "Rates Only" analysis. 
The minimum ^"^^ the corresponding g.o.f. are given in the last two columns. The 
number of degrees of freedom is 3: 5 rates (Homestake, SAGE, Gallex/GNO, SNO, SK) - 2 
parameters. 



Solution 




/eV^ 


tan2 6' 


Amm 


g.o.f. 


LMA 


2.9 


X 


10- 


-5 


0.36 


3.55 


0.31 


SMA 


7.9 


X 


10- 


-6 


1.4 X 10-3 


5.1 


0.16 


LOW 


1.0 


X 


10" 


-7 


0.66 


7.9 


0.05 


VAC 


7.9 


X 


10- 


-11 


3.5 


2.24 


0.52 


Just-so^ 


5.5 


X 


10- 


-12 


2.0 


16.4 


0.0009 


SMA (sterile) 


4.4 


X 


10- 


-6 


1.0 X 10-3 


17.4 


0.0006 


VAC (sterile) 


1.0 


X 


10- 


-10 


0.35 


6.4 


0.094 



2.4 Analysis of "Rates Only" 

To clarify the relative significance of the total rates and the SK spectrum in the global 
analysis of the data we have performed a fit to the rates in the 4 experiments: Homestake, 
SAGE, GALLEX/GNO, SNO and SK. In this analysis we use the total rate of events in SK 
but do not use the SK energy spectrum of the recoil electrons. The results of the test 
are summarized in the Table 4. 

Notice that with the SNO result the SMA solution does not give anymore the best fit 
to the "rates only", in contrast with pre-SNO analyses. The best fit is obtained in the VAC 
solution region. However, parameters of this solution differ significantly from the parameters 
of the global solution when spectral data are included. 

The "rates only" fit shifts the LMA region to smaller Am^. The LOW solution is 
practically unchanged. Now VAC(sterile) and LOW have low goodness of fit. 

Notice that in the absence of the spectral data the value of Xmm is comparable or larger 
than the number of d.o.f. . 

2.5 Pull-ofF diagrams 

In order to check the quality of the fits we have calculated predictions for the available 
observables in the best fit points of the global solutions found in the free fiux analysis (see 
Table 5). Using these predictions we have constructed the "pull-off' diagrams (fig. |^) which 
show deviations, Dk-, of the predicted values of observables K from the central experimental 
values expressed in the la unit: 

Dr = K = Qati Qce, RcCi Rue-, ^^DN- (6) 
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Table 5: Values of the total rates in CI— , Ga—, SK and SNO experiments in the best fit 
points of global solutions found in the free flux analysis. The rates in the radiochemical 
experiments are given in SNU. For SuperKamiokande and SNO the ratios of the best fit 
rates to the rates predicted in the SSM is given. 



Solution 


CI 


Ca 


SK 


SNO 


LMA 


2.89 


71.3 


0.452 


0.323 


SMA 


2.26 


74.4 


0.463 


0.396 


LOW 


3.12 


68.5 


0.446 


0.368 


VAC 


3.13 


70.2 


0.423 


0.364 


Just So^ 


3.00 


70.8 


0.434 


0.434 


SMA(sterile) 


2.93 


75.5 


0.435 


0.445 


VAC (sterile) 
Just So^ 


3.24 


69.9 


0.414 


0.381 


3.01 


70.9 


0.434 


0.435 



Here ax is the one sigma standard deviation for a given observable K. We take the experi- 
mental errors only: ax = cr^^. The theoretical errors are related mainly to the uncertainty 
in the boron neutrino flux. Since is treated as a free parameter we do not take into 
account its theoretical errors. The remaining theoretical errors are small and strongly cor- 
related in Qat, Rcc and R„e- 

The diagrams are good diagnostics of the fit. They allow one to pin down problems that 
some of the solutions have and to elaborate criteria for further checks. 

According to Fig. |^ only the LMA solution does not have strong deviations of predictions 
from the experimental results. LOW, VAC and SMA solutions give somewhat worser fit to 
the data. The fit from other solutions is very bad. 

The pull-off diagrams give some clarification to a common worry, namely that the high 
statistics SK experiment (in particular its spectral data) overwhelms the rates data in the 
global analysis. In particular, solutions, which are strongly disfavored by the rates, give a 
good fit when spectral data are included. 

This problem can be approached in a different way: one can use just one parameter, e.g., 
the first moment, which describes possible distortions of the recoil electron energy spectrum 
25[. 



3 Global solutions: properties and implications 

Here we evaluate the status of global solutions using the following criteria: 
1. Goodness of the global fit in the free flux analysis. 
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2. Deviation of the boron and /lep— neutrino fluxes found in the free flux analysis from 
the SSM values fluxes, that is, the deviation of fs and fhep from 1. 

3. Goodness of the fit in the SSM restricted analysis and in the "Rates only" analysis. 

4. Stability of the solution with respect to variations of the analysis. 

5. Quality of the fit of the individual observables; features of the pull-off diagram. A 
deviation by more than 3a for some observables is a clear signal for trouble. 

We identify solutions which pass all of these criteria. 



3.1 The best fit solution: LMA 



SNO further favors the LMA solution [26 . In all global analyses, in which the SK spectral 



data is included, LMA gives minimal x ■ The best fit point from the free flux analysis is: 

Am^ = 4.8 ■ 10"^ eV^ tan^ 9 = 0.35, fs = 1.2, fhep = 4.0 . (7) 

Large fhep is needed to account for some excess of the SK events in the high energy part of 
the spectrum. For fhep = 1, values of the oscillation parameters are practically the same as 
in the free flux analysis (see Table 2), and the goodness of the fit is even slightly higher. The 
boron neutrino flux is 10% higher than central value in the SSM: Fb = fB'F^^^^ = 5.66-10^ 
cm~^ c~^ being however within la deviation. This flux is rather close to the central value 
extracted from the SNO and SK data 

The fit of the data with the SSM restricted gives minimum of at practically the 
same values of parameters as in (0). 

The values of the oscillation parameters found here are very close to the values found 
by other groups |]|, |^. This shows that the solution is robust and doesn't change with the 
type of analysis. 

Notice that the SNO data lead to a shift of the best fit point (as well as the whole 
region) to larger values of tan^ 6 as was discussed in the introduction. 

According to the pull-off diagram, the LMA solution reproduces observables at ~ la 
or better. The largest deviation is for the Ar— production rate: the solution predicts l.lcr 
larger rate than the Homestake result. 

The "Rates only" analysis shifts the best fit point to smaller Am^. 

Considering the allowed regions at different confidence levels we find the following: 

1) . Am^ is rather sharply restricted from below by the day-night asymmetry of the SK 
event rate: Am^ > 2 • 10"^ eV^ at 99.73% C.L. . 

2) . The upper bound on Am^ is of great importance for future experiments, in particular 



10 



for the neutrino factories |2^. We find from the free flux fit (CHOOZ bound is not included) 



1.9 X 10-^ eV^ 90% C.L. 
Am^ <{ 2.3 X 10-^ eV^ 95% C.L. . (8) 
4.3 X 10-^ eV^ 99% C.L. 



Similar results can be obtained from the analysis in Ref. where also CHOOZ data [28 
have been taken into account. The CHOOZ bound becomes important for larger than 
Am? ~ 8 X 10~^eV^ where it modifies the 3a contour. 

The fit with SSM restricted fs gives stronger bound on Am^. Instead of the limits in 
Eq. (I) we get the upper bounds 1.7 • 10-^ 2.1 • 10-^ and 3.1 • 10"^ eV^ for 90, 95 and 99 
% C.L. correspondingly. 

3) . The SNO and SK results (evidence of z/^, z/^ appearance) give an important lower 
limit on mixing: 

tan2^> 0.2 99% C.L.. (9) 

4) . Maximal mixing is allowed only at the ~ 3a level: 

tan^e> 1 for Am^ = (4 - 10) • 10"^ eV^ 99.73% C.L.. (10) 

In spite of the shift of the best fit point to larger values of Am^ the C.L. for acceptance 
of maximal mixing is not lower than it was before the SNO result. The reason is that the 
SNO rate corresponds to a survival probability smaller than 1/2, which disfavors maximal 
mixing. 

Similar result follows from the analysis in Ref. P where it was found that maximal mix- 
ing (at 3a level) is allowed for Am^ = (4 - 20) • 10"^ eV^. 

In the fit with the SSM restricted Jb, maximal mixing is even more disfavored: tan^ 6 < 
0.9 at 99.73 % C.L.. 



3.2 Large or Small? The fate of SMA 

The fate of the SMA solution, the only solution which is based on small mixing, is of great 
importance for future developments in both theory and experiment. 

We find that the SMA solution does not appear at 3a level (with respect to the global 
minimum) in the free flux fit. The best fit point parameters are: 

Am^ = 6.0 ■ 10"*^ eV^ tan^ 9 = 0.0019, fs = 1.12, fhep = 4. (11) 

SNO shifts the local minimum to substantially larger mixing angles in comparison with the 
pre-SNO result. This is a consequence of appearance of the u^/ut-— flux, which implies large 
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transition probability and therefore large mixing angles. At such a large tan^ 9 one expects 
significant distortion of the boron neutrino spectrum (see sect. 4). 

Important feature of the solution is large fiux of the /lep— neutrinos. It is this large 
fiux which, together with correlated systematic error and > 1, makes possible to get 
a reasonable description of the SK energy spectrum. If the SSM value is taken for the 
/iep— neutrino flux (see Table 2) the increases by Ax^ = 1.4. 

These results are in agreement with those obtained in ^j. In the SMA is accepted 
at lower than 99% level. Surprizingly, in this analysis the SNO result shifts mixing to even 
smaller values: the best fit value from is at tan^ 9 = A - 10~^. 

Before the SNO result, the SMA solution was always giving the best fit to the total 
rates. Inclusion of the CC— event rate measured by SNO moves the SMA solution to third 
position, after VAC and LMA. Furthermore, the fit is no longer good: Xmm = ^-l 3 
d.o.f.. From the pull-off diagram we find that there is a tension between the SNO and 
Homestake rates: The SNO data (CC— rate) requires rather small survival probability for 
boron electron neutrinos. Furthermore Gallium experiments imply strong suppression of 
the Beryllium neutrino flux. This leads to low (1.9(j) Ar— production rate. At the same 
time, according to flg. |^ the CC— event rate is l.Tcr higher than the SNO result. 

For mixing angle corresponding to the best flt point (|ll|) one expects signiflcant regen- 



eration effect in the core-night bin p9|] due to parametric enhancement of oscillations for 
the core crossing trajectories ||3^. However, the day and the night spectra we used in our 
analysis are not sensitive to this feature. The zenith angle distribution of events measured 
by SK does not show any core enhancement Inclusion of this information in the global 
analysis will further disfavor the SMA solution. 



3.3 LOW: next best? 

For the best flt point of the free flux analysis we get: 

Am^ = 1.1 ■ 10"^ eV^ tan^ 9 = 0.66, fn = 0.88, fhep = 2.0. (12) 

If the /lep- neutrino flux is flxed at its SSM value, fhep = 1, the best flt point shifts to larger 
mixing: tan^ 9 = 0.69. Notice that the solution implies ~ la lower boron neutrino flux than 
in SSM. 

In the analysis with SSM restricted fs the best flt point is the same as in ([T^ ) 

The LOW solution gives rather poor flt of the total rates = 7.9 for 3 d.o.f. (see Table 
4). In the best flt point we get 2a larger Ar— production rate and 1.9cr lower Ge— production 
rate. It is this case when the overwhelming spectral information "hides" some problems with 
total rates in the global flt. The use of a single parameter for description of the spectrum 
distortion gives much lower goodness of the flt for the LOW solution as compared with the 
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LMA solution |25| 



3.4 VAC - oscillation solution is back ? 

As anticipated from the comparison of the pre-SNO predictions with the SNO result the 
VAC solution improves its status. Indeed, in the best fit point of the free flux analysis: 

Am^ = 1.4 • 10"%V^ tan^ 6 = 0.40 (2.5), fs = 0.53, fhep = 6.0. (13) 

is even lower than the LOW solution has. The solution with parameters ([T3|) has been 
found in the pre-SNO analysis 0, but before the SNO result the goodness of the fit was 
substantially lower in comparison with other solutions. 

The solution gives very good description of the SK energy spectrum. The is sub- 
stantially smaller than the number of degrees of freedom. The solution reproduces rather 
precisely the bump in both the night and the day spectra at (7 - 8) MeV and the dip at (11 
- 12) MeV. (This features can be well seen in fig. 4e from the Ref. [^.) The bump in the 



spectrum originates from the first maximum of the oscillation probability which corresponds 
to the oscillation phase 27r. Above the bump the probability decreases with energy and the 
increase of R^e at E > 12 MeV is due to large flux of the /lep— neutrinos. 

Notice that an excellent description of the spectrum requires non-maximal mixing, oth- 
erwise the distortion is very strong. Value of sin^ 29 which immediately determines the 
depth of oscillations should be about 0.8. Then to compensate for rather large survival 
probability and to explain the SNO result one needs to assume a small {fs ~ 0.5) original 
boron neutrino flux. 

However, there are several problems with this solution. It requires substantially (~ 3a) 
lower original boron neutrino flux than in SSM and substantially higher original /lep— neutrino 
flux: fhep = 6. In the fit with fhep = 1 the increases by Ax^ = 1.4. 

For this solution one predicts the seasonal asymmetry due to oscillations A ^ —O.GAq, 
where Aq ~ 7% is the asymmetry due to the geometrical factor only (1/-R^ change of the 
flux). Thus one expects a suppressed seasonal asymmetry in contrast with observations. 

The solution gives rather good fit to the rates. Although in the "Rates only" analysis 
the best fit point shifts to a different island in parameter space. According to the pull-off 
diagram, the solution predicts 2. la higher Ar— production rate and 2. 60" lower SK rate, and 
no day-night asymmetry. 

The SSM restricted global fit shifts the best fit point to an "island" centered around 

Am^ = 4.8 • 10-^° eV^ tan^ ^ = 1.9. (14) 

Now fs = 0.72. This solution was found by other groups too. The solution in the same 
"island" of the oscillation parameter space already appeared before (after 508 days of SK 
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operation) when a significant excess of events at tlie liigfi energy part of the spectrum was 
observed. The solution was later excluded by SK data on the spectrum. After the SNO 
result it re-appeared again. The solution can reproduce some bump at ii^ ~ 8 MeV and 
increase of R^e at the high energies. 

Thus, in the VAC region there are three local minima (two of them are degenerate) 
with rather close x^- Small variations in the analysis shift the best fit VAC point from one 
minimum to another. 

Although the VAC solution provides a very good fit to the data in the free flux global 
analysis, it does not pass additional criteria of quality. It requires very strong deviations of 
the boron and hep neutrino fluxes from the SSM values. The goodness of the fit becomes 
substantially worse when the SSM restrictions are imposed on these fluxes. There are sig- 
nificant deviations in the pull-off diagram. 

The Just-so^ solution looks extremely unlikely in the light of the SNO result. It gives 
a very poor fit to the rates. In fact, the fit is so bad that even the fiat spectrum that 
it predicts, in agreement with the SK measurement, is not sufficient to make it plausible. 
Given the excellent fit of the LMA solution, the Just-so^ solution is ruled out at 3a C.L.. 
This result is in agreement with Ref . . Note that if SNO fails to observe a large day-night 
asymmetry, the situation might change and the Just-so^ solution might reappear at about 
~ 3a C.L.. 



3.5 How large is large mixing? 

This question is crucial for theory. In a number of approaches to the bi-large mixing one 
gets mixing of the electron neutrino which is very close to maximal mixing (see |^ for a 
general discussion). 

In the LMA region we find from the free flux analysis: 

r 0.68 95% C.L. 
tan^^ < I 0.82 99% C.L. (15) 
[ 1.05 99.73% C.L. 

and, as we discussed in sect. 3.1, maximal mixing is allowed at ~ 3cr level for Am^ = 
(4 — 10) ■ 10~^ eV^. In the SSM restricted analysis the bounds become stronger. 

In the LOW region we find tan^ < 0.9 at 95% C.L.. Maximal mixing is accepted at 
slightly lower than 99% C.L. . in the range (6 ■ 10~^ — 3 • 10^'') eV^ which covers LOW and 
the so called QVO (Quasi- Vacuum Oscillation) regions. 

Maximal mixing is the best fit value of the Just-so^ solution. Although this solution is 
ruled out at 3a level in the global analysis. Notice, the fit of the data in the VAC region 
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implies significant deviation from maximal mixing. 



3.6 Do pure sterile solutions exist? 

The best (and the only accepted at 3a level) pure sterile solutions is the VAC (sterile), as 
could be expected from our pre-SNO analysis [Q. In the free flux analysis the solution 
appears at 90% C.L. already with the best fit point: 

Am^ = 1.4 ■ 10"^° eV^ tan^ 9 = 0.38 (2.6), fs = 0.54, f^ep = 9. (16) 

Partially, the difference between the SK and SNO rates is explained by distortion of the 
spectrum: the suppression increases with energy and therefore the higher threshold at SNO 
leads to lower averaged survival probability. However, mainly the fit "shares" the deviations 
between SNO and SK: the solution requires 1.2a higher SNO rate and significantly lower 
boron neutrino flux which gives a ue— scattering rate 3. Oct below the SK result (see fig. 
The free flux analysis without SK total rate does not locate this problem immediately. In 
principle, it should reappear in the fit of the spectrum with similar Ax^ since the solution 
fits the shape of the spectrum rather well. 

Notice that the parameters of this solution are very close to the parameters of VAC (active) 
apart from the two times larger hep— fiux. So this solution leads to the same type of the 
spectrum distortion as the one described in sect. 3.4 with the bump at i^^ = 7 — 8 MeV and 
the dip at E = 11 - 12 MeV. 

In addition to a strong deviation of the ue scattering event rate from the SK result, the 
solution predicts 2. la larger Ar— production rate. The fit worsens when SSM restrictions are 
imposed. For the SSM value of the /lep— neutrino flux the best fit point shifts to a smaller 
mixing angle and the increases by Ax^ = 1.8. In the analysis with SSM restricted Jb, 
the goodness of the fit drops further: Ax^ = 6.2. In the fit of the rates the solution shifts to 
smaller Am^ where a description of the SK spectra becomes bad. Thus, the solution does 
not pass additional quahty tests. 

The SMA(sterile) solution gives a very bad fit since it leads (in contrast with the VAC 
solution) to a distortion of the boron neutrino spectrum with suppression which weakens 
with increase of energy. This solution contradicts the SNO result: 3.3a higher rate is pre- 
dicted. 



Concerning sterile solutions one remark is in order. Better fit of the data can be obtained 
if more than 1 sterile neutrino participate in the conversion. Such a possibility is realized 
when solar neutrinos convert to the so called "bulk" neutrinos which propagate both in 
usual and in (large) extra space dimensions [Q. From the four dimensions point of view 



the solar is transformed to several Kaluza-Klein mode of the bulk neutrino which show 
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up as sterile neutrinos. In this case one can reproduce (with low x^) the SK spectral data, 
and at the same time have better descriptions of the total rates (see for recent analysis). 



4 SNO: Predictions for the next step 



Forthcoming results from SNO will include measurements of the day-night asymmetry and 
a more precise determination of the electron energy spectrum (higher statistics and lower 
energy threshold). Later results on the NC/CC ratio will be available. Predictions for these 
observables have extensively been discussed before [^5|, ^ |39|, ^ . Here we sharpen 
the predictions using the latest solar neutrino data. We also calculate the expected values 
of the total event rate in the BOREXINO experiment which start to operate soon. 



4.1 Day-Night asymmetry 

We calculate the D-N asymmetry at SNO, Afy^, defined as 

^DN -^N + D ^^^^ 

for the events above the threshold E^^ = 6.75 MeV. 

We compare the asymmetry of the CC— events at SNO with the asymmetry of the 
i/e— events measured at SK: A^^. There are tree factors which can lead to substantially 
different SNO and SK asymmetries: 

A "dumping" factor, rj^ump, describes the suppression of the D-N asymmetry in the z^e— 
event rate due to the contribution of the z/^, scattering to the signal. We get |T6[ 



■^DN ^ VdumpAjyj^, (18) 

where 

Vdump = 1 + -^Y^Typ • (^9) 

Here r = cr(z/^)/(T(z/e) is the ratio of cross-sections of the t'^e— and t'eC— scattering, and 
P is the averaged survival probability. Notice that with decrease of P the damping factor 
increases. 

A second factor, rjthr, describes the effect of difference of the energy thresholds: 5 MeV 
for SK and 6.75 MeV for SNO. It also accounts for the larger minimum difference (the 
binding energy of the deutron, 1.44 MeV) between neutrino energy and electron energy in 
SNO. 
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A third factor, rjreg, is related to the difference of the geographical latitudes. Since the 
Sudbury mine is at higher latitude, the regeneration effect there is slightly weaker than at 
Kamioka. 

The total difference between the SNO and SK asymmetries is the product of these three 
factors. 



Let us consider now the predictions for the asymmetries for individual solutions. 

1 ). LMA- solution. In Fig. ^ we show dependence of the asymmetry on Am^ for different 
values of tan^ 6 from the allowed region. The asymmetry decreases with increase of Am^ as 
~ 1/Am^ for Am^ > 4 • 10~^ eV^ and at larger Am^ the decrease is faster due to effect of 
the adiabatic edge (see for details). The dependence of asymmetry on tan^6' is rather 



weak. In the best fit point we get 

Af)%^ = 7.2% {E'^ = 6.75 MeV). (20) 

In the allowed region the asymmetry can take any value from practically zero to 15% at 90 
% C.L.. At 3a level it can reach 20 %. The asymmetry is maximal at the smallest possible 
values of Am^ and tan^ 6 within the allowed region. It increases with the energy threshold 
due to an increase of the regeneration factor: freg oc E/Am"^. 

The expected asymmetry at SNO is substantially larger than at SK. In the best fit point 
we get Af)^ = 3.6%, thus A^^ ^ 2 • Af,^. This difference can be easily understood con- 
sidering the factors mentioned above. Indeed, for the best fit point the survival probability 
can be estimated as P ~ Rsno/Jb- Inserting numbers from the Table 5 we get from Eq. 
([T9| ) rjdump = 1-6. In the LMA region the asymmetry increases with Eth- The threshold 
factor equals rithr ~ 1.2, and rjreg is close to 1. As a result, the overall enhancement factor 
at SNO is about 2. 

2). LOW solution. The dependence of the asymmetry on Am^ for different values of 
tan^ 9 from the allowed region is shown in Fig. ^. In the best fit point (of the free flux 
analysis) 

Ag)^^ = 2.4% {E'^ = 6.75 MeV). (21) 

The asymmetry increases linearly with Am^ for Am^ > 10"'' eV^. It can reach 10 - 12 % at 
the upper border of the allowed {3a) region. For Am^ < 10"'' eV^ the asymmetry decreases 
faster with Am^ due to effect of the non-adiabatic edge (see fl^ for details). It depends 
weakly on the mixing angle. 

Let us emphasize that the D-N asymmetry in the best fit point of the LOW solution is 
substantially smaller than the one in the LMA region. Thus, an observation of asymmetry 
AfJ^ > 5% will strongly favor the LMA solution. 

In the LOW region the D-N asymmetry at SNO is also larger than in the SK detector. 
However the difference here is not as large as in the LMA case. In the best fit point we get 
= 2%, which results in Af)%^ = l.2A^J^. There are two reasons for such a difference: 
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(i) The average survival probability is larger now (mixing angle is larger) : P ^ Rsno / fs = 
0.42 (see Table 5). As a consequence, the damping factor is smaller: rj^amp = 1-45. 

(ii) The difference in the thresholds works in the opposite direction in comparison with 
the LMA case, thus suppressing the SNO asymmetry. Indeed, in the LOW region the regen- 
eration factor decreases with E/Am'^ and therefore the asymmetry decreases with increasing 
energy threshold. For rjthr ~ 0.85 we get a total enhancement factor in agreement with the 
exact calculation. 

3). SMA solution. The dependence of the asymmetry on tan^^ for different values 
of Am^ is shown in Fig. ^ In the best fit point we get Afj%° = 2.6% and most of this 
asymmetry is collected from the "core" bin. The asymmetry increases fast with tan^ 6. In 
the 3a allowed region it can be as large as 8%. 

The SNO asymmetry is only slightly higher than the SK- asymmetry: = 2.4%. In 

the SMA range due to strong dependence of the survival probability on energy the relation 
between the SNO and SK asymmetries is more complicated than in Eq. ([T9|). 



4.2 Spectrum distortion 



As follows from previous studies [^, in general the sensitivity of the SNO measure- 

ments to the energy spectrum distortion is higher than the sensitivity of SuperKamiokande 
due to better correlation of the neutrino and the (produced) electron energies. However, 
the present statistics in SNO is much lower than in SK, and the statistical errors are more 
than 2 times larger (in the range 8-10 MeV we get 10 -12 % at SNO as compared with 4 - 
5 % at SK). 

In our analysis we do not use the spectral information from SNO. Instead, we present 
here a qualitative discussion comparing the present SNO data with predicted spectra from 
different solutions. This allows us to evaluate significance of the present and forthcoming 
SNO results. 

In Fig. ^ we show the expected spectra of events at SNO for several values of the 
oscillation parameters in the LMA region. In the high energy part of the spectrum the 
distortion is due to the earth regeneration effect as well as contribution of the hep— neutrino 
flux. However, for Am^ < (4 — 5) • 10~^ eV^ the regeneration effect is small and the turn 
up at > 11 MeV is due to increased {fhep = 5) flux of the /lep— neutrinos. The turn 
up of the spectrum at low energies is the effect of the adiabatic edge of the suppression 
pit. Recall that the ratio of the event rate with oscillations to the event rate with no 
oscillations is determined by the product of the survival probability times the factor fs'- 
Rcc oc fB{A.rn?)P[Aim?). With increase of Aw? spectrum shifts to the adiabatic edge; for 
Air? < 10^^ eV^ the probability can be written as 
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where A is the correction due to adiabatic edge and it is proportional to the first moment 
5T: 

Aoc5Toccos2^ (22) 



(see Appendix of [|T^ for more details). The factor cos 26' accounts for the disappearance 
of the distortion when the mixing approaches maximal value. With increasing Am? the 
size of the turn up and the distortion of the spectrum (first moment) increase as (Am^)^. 
Notice that at the same time decreases to compensate for a total increase of the survival 
probability. The distortion reaches a maximum at Am^ ~ 1.5 • 10~^ eV^ when the boron 
neutrino spectrum is in the middle of the adiabatic edge. With further increase of Am? , the 
spectrum shifts out of the suppression pit and the distortion decreases. For Am^ > 3 • 10~^ 
eV^ the spectrum is in the range where the probability is determined basically by the 
averaged vacuum oscillations: P = 1 — sin^26'/2 and does not depend on the neutrino 
energy. Notice that the distortion of the spectrum weakens due to integration over neutrino 
energy and folding with the energy resolution function. Moreover, the SNO sensitivity to 
the distortion at low energies is weakened by the fast decrease of the cross-section with 
energy. 

The curves shown in Fig. ^ illustrate this behavior of the spectrum. The short dashed 
line corresponds to near maximal distortion. As follows from the figure it will be very 
difficult to establish the distortion expected from the LMA solution with SNO. One needs 
to measure the spectrum down to ~ 5 MeV and to increase substantially the statistics. 
Thus, the prediction for SNO is that no significant distortion should be seen in forthcoming 
measurements. 

In Fig. 1^ we show the expected spectra of events at SNO in the best fit points of the 
LOW, SMA and VAC (active) solutions. The correlated systematic errors (mainly due to 
the error in the absolute energy scale calibration) are not shown here. These errors can 
affect the conclusions from the fit to the spectrum, making the spectrum appear fiatter. 

In the best fit point of the LOW solution the Earth regeneration effect on the shape of 
the spectrum is small. The weak positive slope (positive shift of the first moment) is due 
to the effect of the adiabaticity violation (non-adiabatic edge of the suppression pit). The 
slope increases with decreasing tan^ 6 as well as Am^ (for some analytical studies see JT^)- 
It will be difficult to establish this distortion with SNO. 



For the VAC (active) solution with low Am^, Eq. (p!3D, one predicts maximum of the 
ratio Rcc at ~ 6 MeV. It corresponds to the first oscillation maximum of the survival 
probability. The ratio decreases with increase of energy and at ii^ > 11 MeV the distortion 
fiattens due to contribution of the /iej9— neutrino fiux. Thus, negative shift of the first 
moment (slope) is expected. As follows from the figure, the distortion is at the border of 
sensitivity of already present SNO data (correlated systematic errors can partly improve 
an agreement between the data and the predictions). Further increase of statistics and 
especially measurements of the spectrum at lower energies will be crucial for discrimination 
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of the solution. Notice that the oscillation parameters are well fixed by the SK spectrum, 
and no freedom exists, e.g., to change the distortion by varying Am^. In particular, the 
position of the maximum at ii^ ~ 6 MeV is immediately determined by the maximum in 
the SK spectrum at 8 MeV. The distortion can be identified by comparison of the averaged 
ratio Rcc at low and high energies: Rcc{< 9 MeV) and Rcci> 9 MeV). 

Similar distortion is expected for the VAC(sterile) solution, see Fig. 7c. For VAC (active) 
with high Am^ (see Eq. (H)) one expects fiat distribution at low energies due to strong 
averaging effect and a bump at high energies E > 12 MeV. 

For the SMA (active) solution there is a strong positive shift of the first moment 
(slope). Correlated systematic errors improve the agreement with the experimental data. 
The present SNO data can already have some impact on this solution. The allowed region 
will drift to smaller mixing angles where the distortion is weaker. The best fit boron fiiix 
should decrease in order to compensate for the increase of the survival probability. 

For SMA(sterile) the distortion is very weak, as a consequence of very small mixing angle. 
However for this solution the predicted spectrum is systematically above the experimental 
points which corresponds to 3.3a higher total rate. 

In conclusion, the present SNO spectral data further favor the LMA and LOW solutions. 
However, it will be difficult to establish with SNO the weak spectral distortions expected 
for these solutions. At the same time the SNO measurements can be sensitive to strong 
distortion of the spectrum predicted by already disfavored solutions such as VAC and SMA. 



4.3 NC/CC ratio 

The reduced neutral current event rate [NC] is defined as the ratio of the rates with and 
without oscillations: [NC] = N^c / N^q^ . Similarly, the reduced rate of the charged current 
event rate equals [CC] = Ncc /N^^^^ . We have calculated the ratio of the NC and CC 
reduced rates, [NC]/[CC], for global solutions found from the free fiux fit. 



For the active neutrino conversion the ratio equals 

NC 1 fB 



CC P R 



SNO 



(23) 



where P is the effective (averaged) survival probability of the electron neutrinos. In the 
best fit points, using results for fs and Rsno from the Tables 1 and 5, we find [NC]/[CC] 
= 3.5 (LMA), 2.4 (LOW), 2.8 (SMA), L5 (VAC), ^ 1 (Just-so^) in agreement with resuhs 
of numerical calculations (see fig. H). 

For the sterile solutions we have: 

[CC]- Pec' ^ ^ 
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where Pnc and Pec are the effective survival probabihties for the NC and CC samples 
correspondingly. For the NC sample the energy threshold is about 2.2 MeV, so that P^c is 
averaged over larger interval than Pec- As a consequence, the probabilities Pec and Pnc 
are in general different. 

Notice that the largest value of the ratio is expected for the LMA solution: 

M = 3.5l- (LMA). (25) 

Then rather close predictions follow from the SMA and LOW solutions. Smaller value is 
predicted for VAC. For sterile solutions the ratio is close to 1. 

According to fig. ^, there is a significant overlap of the 3cr regions of predictions. How- 
ever, measurements of the ratio with better than 20 % accuracy can significantly contribute 
to discrimination of the solutions. 



4.4 BOREXINO rate 

For global solutions we have calculated the reduced event rate in the BOREXINO experi- 
ment Rborexino = Nosc/^ssM, whcrc Nqsc and NssM are the expected rates with and 
without neutrino conversion. We averaged the effect over the year. 

The rate can be estimated as Rborexino = Psei^ — r) + r, where P^e is the averaged 
survival probability for the ^Be- neutrino line. The probability Pbs has been calculated for 
the oscillation parameters from the 3a allowed regions of solutions found from the free flux 
analysis, r is the ratio of cross-sections of the z/^e— and z/ge— scattering at the the beryllium 
line. We have used the neutrino-electron scattering cross-sections from Ref. . 



The results of calculations for six global solutions are shown in the Fig. The rate is 
suppressed for the solutions. In particular, for the LMA solution we get 

Rborexino = 0.66^o;Qg (LMA). 

Similar suppression is expected for the LOW. The VAC solutions predict stronger suppres- 
sion, the lowest rate is expected for SMA. Notice that the range of predictions for SMA is 
very small due to smallness of the allowed region itself. 

There is a significant overlap of the predicted intervars for LMA, LOW and VAC. There- 
fore, it will be difficult to discriminate among these solutions using just total rate measure- 
ments. However, Borexino has other capabilities, e.g. it can detect strong day-night effect 
in the case of the LOW solutions whereas strong seasonal variations are expected if 



QVO or VAC solution is realized |45 
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5 Conclusions 



We have performed a global analysis of the solar neutrino data including the charged current 
event rate measured by SNO. We tested the robustness of the global solutions by modifying 
the analysis. The quality of the fit was checked by construction of the pull-off diagrams. 

We find that the LMA solution with parameters in the range Am? ~ (4 — 5) • 10^^ eV^ 
and tan^ 9 — 0.35 — 0.40 gives the best fit to the data. Moreover, the solution reproduces 
the flat spectrum of recoil electrons in SK and gives very good description of the rates and 
of the Day-Night asymmetry. It is in a very good agreement with SSM fluxes of the boron 
[fs = 1-10 — 1.13) and /lep— neutrinos. The vahics of the oscillation parameters and the 
goodness of the fit arc stable with respect to variations of the analysis. 

The LOW solution appears at 90 % C.L. with respect to the best fit point in the LMA 
range. It gives good fit to the SK spectral data, but rather poor fit to the rates: it predicts 
larger than measured 74r— production rate and smaller than measured Ge— production rate. 

The VAC solution has high goodness of the global fit due to a very good description 
of the day and the night spectra measured by SK. It has, however, problems with other 
criteria. The solution requires strong deviation of the boron and /lep— fluxes from their 
SSM values. The fit becomes substantially worse when SSM restrictions arc applied. The 
solution shows strong deviations in the pull-off diagram (especially for the SuperKamiokande 
rate and the 74r— production rate). The analysis of the rates only selects a different point 
in the oscillation parameter space. 

The SMA appears at about 3(T C.L. level with respect to the best global solution in the 
LMA range. It gives bad fit of the recoil electron energy spectrum and also the fit of the 
rates is rather poor. 

The best sterile solution is VAC(sterile) with Am^ = 1.4-10"^*' cV^. It is the only sterile 
solution accepted at lower than 3(7 C.L. in the free flux analysis. This solution, however, 
has serious problems with other tests. Similarly to VAC (active) the solution requires strong 
deviation of the boron and hep— fluxes from their SSM values. The fit becomes substantially 
worse when SSM restrictions are apphed. The solution shows strong deviations in the 
pull-off diagram: in particular, the ve— event rate is about 3o" below the SK rate, the 
Ar— production rate is 2Aa higher than the Homestake result. The analysis of the rate 
only selects a different point in the oscillation parameter space. Significant distortion of the 
energy spectrum is expected. Appearance of the VAC solutions (both active and sterile) in 
the global fit is related to a small spread of the experimental points in the SK spectrum. 
These solutions can describe rather precisely the bump at 7 - 8 MeV and the dip at 11 - 12 
MeV. 

Maximal mixing is allowed at 3o" level in the LMA region and at 99 % C.L. in the LOW 
and quasivacuum oscillation regions. 

Measurements of the D-N asymmetry will provide strong discrimination among the 
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solutions. In the best fit range of the LMA solution the asymmetry in SNO is ~ 
(7 — 8)%, and in the Sex allowed region it can reach 15 - 20 %. In the LOW region one 
expects smaller asymmetry: in the best fit point ~ (2 — 3)%, although in the 3(T 

allowed region it can be as large as 10 - 12 %. For the SMA solution the asymmetry is 
Af)%^ < 3%. The predicted SMA zenith angle distribution is not supported by SK data. 

Clearly, observation of a large D-N asymmetry A^^ > 5% will favor the LMA solu- 
tion. It will strongly disfavor the SMA solution, and exclude solutions based on vacuum 
oscillations (VAC, QVO and and Just-so^). Furthermore, this result will strongly restrict 
LMA to lower Am^. Even with existing statistics (expected error is 3 - 4 %) any SNO result 
on the D-N asymmetry will be of physical importance excluding some solutions or further 
restricting the allowed regions. 

Significant distortion of the boron neutrino spectrum is expected for SMA and VAC 
solutions and already the present data can affect them. The VAC solution predicts the 
bump at 6 MeV and the dip at 11 MeV in the dependence of the Rcc on energy. In the case 
of LMA solution turn up of the spectrum at low energies is expected which will be rather 
difficult to observe with SNO. For LMA and LOW one expects practically no distortion in 
future SNO measurements. 

Important discrimination among solutions can be done using precise determination of 
the [NC]/[CC] ratio. 

With forthcoming SNO data, we have a chance to make further major step in identifi- 
cation of the solution of the solar neutrino problem. 
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Figure 1: Global solutions from free flux analysis. The boron and hep— neutrino fluxes are 
considered as free parameters. The best fit points are marked by dark circles. The absolute 
minimum of the is in the LMA region. The allowed regions are shown at 90%, 95%, 99% 
and 99.73% C.L. with respect to the global minimum in the LMA region. 
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Figure 2: Global solutions for the SSM restricted boron neutrino flux, hep— neutrino flux 
is free. The best fit points are marked by dark circles. The absolute minimum of the is 
in the LMA region. The allowed regions are shown at 90%, 95 %, 99% and 99.73% with 
respect to the global minimum in the LMA region. 
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Figure 3: Pull-off diagrams for global solutions. Shown are deviations of predictions for 
the y4r— production rate Ge— production rate, SK rate, the day-night asymmetry at SK, 
and the SNO rate from experimentally measured values. The pull-offs are expressed in the 
units of 1 standard deviation, la. 
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Figure 4: The dependence of the Day-Night asymmetry of the CC— event rate measured 
at SNO on the Am^ for different values of the mixing angle. The oscillation parameters are 
taken from the LMA allowed region. The best fit value of Am^ is marked by the vertical 
dotted line. 
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Figure 5: The same as in the Fig. 3a for the LOW solution. The dotted vertical lines mark 
the best fit values of Am^ from the free flux fit (left) and SSM restricted fit (right). 
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Figure 6: The dependence of the Day-Night asymmetry of the CC— event rate measured 
at SNO on tan^ 6* for different values of Am^. The oscillation parameters are taken from 
the SMA allowed region. The best fit value of tan^ 6 is marked by the vertical dotted line. 
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Figure 7: The recoil electron energy spectra of the CC— events in SNO for global solutions 
of the solar neutrino problem. Shown are the ratio of the number of events with and without 
conversion as a function of the electron kinetic energy, a). The spectra in the LMA solution 
region for different values of Am^ and tan^6' = 0.35. b). Spectra for the best fit points 
of the SMA, LOW and VAC (active) solutions, c). Spectra for the best fit points of the 
SMA(sterile) and VAC(sterile) solutions. Shown also in all panels is the SNO experimental 
data. 
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Figure 8: The ratio of the reduced NC - event rate to the reduced CC- event rate. The 
circles give values of [NC]/[CC] in the best fit points of the global solutions. The error bars 
show the prediction intervals of [NC]/[CC] which correspond to the 3cr allowed regions of 
global solutions. 
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Figure 9: The reduced event rate in the BOREXINO experiment. The circles give values 
of Rborexino in the best fit points of the global solutions. The error bars show the prediction 
intervals which correspond to the 3a allowed regions of global solutions. 



35 



